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ABSTRACT 

We present the morpho logical catalog of galaxies in nearby clusters of the WINGS 
survey ( [Fasano et al. 2006). The catalog contains a total number of 39923 galaxies, 
for which we provide the automatic estimates of the morphological type applying the 
purposely devised tool MORPHOT to the V-band WINGS imaging. For -3000 galax- 
ies we also provide visual estimates of the morphological types. A substantial part of 
the paper is devoted to the description of the MORPHOT tool, whose application is 
limited, at least for the moment, to the WINGS imaging only. The approach of the tool 
to the automation of morphological classification is a non parametric and fully empiri- 
cal one. In particular, MORPHOT exploits 21 morphological diagnostics, directly and 
easily computable from the galaxy image, to provide two independent classifications: 
one based on a Maximum Likelihood (ML), semi-analytical technique, the other one 
on a Neural Network (NN) machine. A suitably selected sample of —1000 visually clas- 
sified WINGS galaxies is used to calibrate the diagnostics for the ML estimator and as 
a training set in the NN machine. The final morphological estimator combines the two 
techniques and proves to be effective both when applied to an additional test sample 
of —1000 visually classified WINGS ga laxies and when compared with small samples 
of SDSS galaxies visually classified by 'Fukugita et al.| ( [2007 1) and ^ir et al.| (|2010j). 



Finally, besides the galaxy morphology distribution (corrected for field contamination) 
in the WINGS clusters, we present the ellipticity (e), color (B-V) and Sersic index (n) 
distributions for different morphological types, as well as the morphological fractions 
as a function of the clustercentric distance (in units of i^20o)- 

Key words: galaxies: clusters - galaxies:general - galaxies: elliptical and lenticular 
- galaxies :cD 



1 INTRODUCTION 

The Wide-field Nearby Galaxy -clusters Survey (WINGS, 
Fasano et al. 2006) has gathered wide-field, photometric 



data, in the optical bands B and V (Varela et al. 2009), 



of several hundred thousand galaxies in the fields of 76 
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nearby clusters (0.04^z^0.07), selected from three X-ray 
flux limited samples compiled from ROSAT All-Sky Survey 
data dEbeling et al.||1996[ [19981 12000| . The observations in 
the optical bands have been obtained with the WFC@INT 
and with the WFI@MPG/ESO-2.2 cameras for the northern 
and southern clusters, respectively. Follow-ups of the opti- 
cal WINGS survey include medium-resolution, multi-fiber 
spectra of ^-^6000 galaxies in 48 WINGS clusters ( Cava et al 
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2009 ), wide-field imaging in the NIR bands J and K of ^10' 
galaxies in 28 WINGS clusters ( [Valentinuzzi et al.|2009 ) and 



U-band, wide-field imaging of 18 WINGS clusters (Omizzolo 
et al.|2011t . Lastly, a narrow H^-band photometric survey is 
presently ongoing on a sub- set of the WINGS cluster sample. 

The WINGS survey was conceived in 2000, mainly with 
the aim of making up for the odd situation for which, while a 
large amount of high quality morphological data for distant 



clusters were already available from HST imaging (,Couch 


et aL||1994| |Pascarelle et aL| 


19951 lOemler et al.||1997| |Kel- 


son et aL||1997 Couch et al. 


1998| Lubin et al.||1998"), high 



quality CCD data were almost lacking for large samples of 
nearby clusters. Actually, the very selection of the WINGS 
cluster sample, as well as the choice of the telescopes and 
the observational constraints of the optical survey, were per- 
formed in order to meet the requirements needed by the 
main original task (morphology of galaxies in clusters), in 
terms of absolute field of view (>1.6 Mpc) and spatial res- 
olution (1'' <1.3Kpc). 

A recent, comprehensive review of the various aspects 
and issues linked to galaxy morphology can be found in B uta| 
(2011). Till a dozen years ago the morphological type esti- 
mate of galaxies has been obtained just by visual inspec- 
tion of photografic plates or CCD frames. A few attempts 
were actually made in the 90's to obtain automated morpho- 
logical classification of galaxies (with neural networks and 
self-organizing maps; Naim et al.||1995 1997| , but they re- 
mained isolated. In the last decade, the sudden availability of 
CCD mosaics has made it no longer feasible to conduct mor- 
phological classifications by eye, since one has often to deal 
with wide and/or deep fields, each one containing thousands 
galaxies. This has triggered a number of papers proposing 
different tools for automatic morphology estimate of large 
galaxy samples. 

There are basically two alternative approaches to the 
problem of automatic morphological classification of galax- 
ies. The first one exploits the parametrization of their radial 
light profiles (see for instance Gutierrez et al.|[2004 Tru- 
illo and Aguerri 2004; Saintonge et al. 2005; Tasca and 
White 2005; Orndahl and Ronnback 2005; Ravindranath et 



al. |2006| [Cassata et al.| 2010). In this case, the most com- 
monly used morphological diagnostics are the bulge frac- 
tion (B/T) and the Sersic's index (n). Several tools have 
been devised to obtain, in (semi) automatic mode, the best- 
fit parameters of the analytical laws used to represent the 
light distribution of galaxies. Among the others, we mention 
GIM2D (|Simard||1998|), GALF IT ( jPeng et al.|[^0Q2t, GAS- 
PHOT dPTgnatelh et aLl|2006|), 2DPH0T (|La Barbera et 



aL] |2008|), GASP2D ([Mendez-Abreu et al.||2008| ) and GAL- 



PATH (Yoon et al. 2009| . However, this approach to the 



morphological classification of galaxies presents two serious 
drawbacks: (i) the analytical components derived from the 
formal best-fitting of galaxy light profiles (usually exponen- 
tial and Sersic laws) often do not correspond to real physi- 
cal components (disk and bulge; see for instance Tasca and 
|White|2005| ); (ii) the correlations between these diagnostics 
{B/T and n) and the visual morphological type are weak 
and show a high degree of degeneracy, especially for early- 
type galaxies (see Figure Vu see also Sanchez- Portal et al, 



)06| ) 



2004 Pignatelh et al. 2006}. These drawbacks reflect the 



fact that structure and morphology of galaxies are intrin- 
sically different concepts (see van der Wei 2008). In fact. 
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Figure 1. Visual morphology versus Sersic's index (top panel) 
and bulge fraction (bottom panel) for the 527 galaxies in common 
between the MORPHOT calibration sample (see Section|3.1J) and 
the WINGS-GASP2D sample ( ,Sanchez-Janssen et al.,201i; . 



while the first one is a global property that can be described 
by means of simple analytical laws and leaves mostly aside 
the problems connected with image texture, S/N ratio and 
resolution (apart from the convolution with the local PSFs), 
the second one mainly deals with pixel-scale behaviours and 
features which, in default of human eyes, make difficult any 
quantitative description. 

The alternative, non-parametric approach tries to face 
the problem relying upon various diagnostics, directly com- 
putable from the digital post-stamp images of galaxies, 
which are empirically found to correlate with the visual 
morphological estimates. The non-parametric tools can be 
in turn divided in two main categories: (i) those using a 
few (two or three) diagnostics and the relative two- or 
three-dimension space to try to segregate galaxies with dif- 
ferent morphological types; (ii) those using Neural Net- 
works (or some other sharp methodology) to combine 
many diagnostics, thus drawing a final, quantitative esti- 
mate of the morphological type. Among the tools belong- 
ing to the first category, besides the pioneering diag- 
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nostic devised by Abraham et al. ( 1996[ concen 



tration vs. asymmetry) and the popular CAS (Concen- 



tration/ Asymmetry /clumpinesS) diagnostic set (Conselice 



2003), it is worth mentioning those proposed by I Abraham 



et al . (2QQ 3t Gini Coefficient ) and |Lotz et"ar] ( |2004[ M20 



coefficient), Lauger et al. (2005 concentration and asymme- 
try at different wavelenghts) , Yamauchi et al.|(|2005| co ncen- 
tration and coarseness coefficients), Menanteau (2006|), 



der Well (|2008) and Petty (2009). To the second c ategory 



belong the tools devised by Odewahn et al.|(2002| Fourier 



(2004 



(2007 



analysis of the images), Ball et al. ( 2004 )7[Coderya et al 



shapelet analysis), Moore et al. (2006), Scarlata et al. 



( 20041 ), |de~la Calleja and Fuentes] (| 20041, [Kelly and McKay 



ZEST), [Huertas-Company et al.| ( ,2008) and [Shamir 



(2009). A mixed approach (5/T decomposition + non para- 



metric diagnostics) has been tried by Rahman and Shan- 
[darin[( [2004t , [Cheng et al.[ ( [2011[ ) and [Vikram et al.[ ( |2010" 
PyMorph). 

The non-parametric approach seems to be more effec- 
tive than the parametric one in estimating the morphological 
type of galaxies ( Hatziminaoglou et al.[[2005 ) and has been 
claimed to be even able to compare with visual estimates 
as far as the intrinsic scatter and the robustness of the re- 
sults are concerned ([Odew ahn et al.[ [2002[ [Beh et al.[[2004 



Huert as- Company et al. [2008). However, a common limita- 
tion of the non-parametric tools available in the literature 
is the scarce ability of separating SOs from elliptical galax- 
ies, which is actually an important issue when dealing with 
galaxy evolution in clusters (iDressler et al."1997! "Fasano et 



aL][20 00; Treu et al. 2003; Postman et al. 2005; Desai et al. 
2007[ [Poggianti et al.[|2009| ). 

In this paper we describe a new automatic, non- 
parametric tool for the morphological type estimate of large 
galaxy samples (MORPHOT Fasa no et "aLl[2007[ ) , which is 
in fact able to separate Es from SOs in a majority of cases. 
MORPHOT belongs to the second previously mentioned 
category of non parametric tools. It starts with a set of 21 
suitably devised morphological diagnostics, and combines 
them in two different (independent) ways, thus producing 
the final morphological type (and the relative confidence in- 
terval) for each galaxy in a given input catalog. We fine-tune 
MORPHOT for extensive application to the WINGS cluster 
sample and present the catalogs of the survey, which contain 
morphological types of ^-^40000 galaxies. We stress that, al- 
though the basic methodology is robust for any set of digital 
images of similar spatial resolution and dynamic range, at 
this stage the tool does not pretend to have a general valid- 
ity, regardless of the observing conditions (telescope, detec- 
tor, seeing) and the galaxy sample (redshift). However, we 
will show that it produces reliable results for the particular 
purposes of the WINGS survey. 

In Section [2] we report about the intrinsic reliabity of 
the visual morphological classifications. In Section |3] we de- 
scribe in some detail the structure of MORPHOT and the 
various steps of the tool's flow-chart. In Section^ we anal- 
yse the performances of MORPHOT on the WINGS galaxy 
sample. In Section |5] we apply MORPHOT to the WINGS 
cluster galaxies, present the WINGS catalogs of morpholog- 
ical types and briefly discuss the main statistical proper- 
ties of galaxy morphology in the WINGS clusters. Section [6] 
summarizes the results and outlines the future employment 
of the MORPHOT classifications. Throughout the paper 



we adopt the following cosmology: Hq—70 kms ^Mpc ^, 
^^=0.3 and Qa=0.7. 



HOW RELIABLE IS THE VISUAL 
CLASSIFICATION ? 



After the pioneering attempt by Raynolds (1920) 



to provide a morphological classification of spi- 
ral nebulae and since the first, definite understand- 
ing by Edwin 'Hubblej ( 1925 ) of the extragalactic na- 
ture of many nebulae, a number of different classifica- 
tion schemes have been proposed for the galaxy morphol- 
ogy. The original Hubble's sequence (Hubble 1922 spi- 



rals/elongated/globular/irregulars) was improved by the au- 
thor himself, first introducing the concept of tuning fork to 
distinguish between normal and barred disk galaxies (Hub- 



ble 



ble 



1926), then defining the SO morphological type (Hub- 



1936 



Later on, the Hubble system has been refined 
and completed, introducing spiral types later than Sc (Sd 
and Sm; de Vaucouleurs 1959), a new type of amorphous 



galaxies (Irr-II; Holmbergyl950|) and ring-based (jde Vau- 



couleurs 1963) or arm-based (Elmegreen et al. 1987) dis- 



tinctions among disk galaxies. 

Radically different classification schemes were proposed 
by [Morgan! ( |1958| and |van den Bergh| ( [19591 |1960a|b| . The 
first one links morphology with central concentration of light 
and stellar populations, also introducing the new cD type. 
The second one links morphology with total luminosity (lu- 
minosity classes) and extends the basic Hubble's scheme 
(E/Sp/Irr) to the lowest luminosity galaxies (dwarfs). 

Today the most frequently used classification scheme 
for statistical studies is the numerical code associated to the 
so called Revised Hubble Type {Trh hereafter), first intro- 



duced by ^de Vaucouleurs | ( |1974| , subsequently improved in 



the Third Reference Catalog of Bright Galaxies (de Vau- 
couleurs et al. 199H RC3 hereafter) and schematically re- 



minded in Columns 1 and 2 of Table [T] 

For reasons which will become clear in the next Sec- 
tion, the MORPHOT tool uses a slightly modified version 
of the Trh code (reported in Column 3 of Table [l] as Tm- 
MORPHOT Type), in which the code -6 is associated to the 
cD galaxies (rather than to compact Es) and the transition 
class E/SO is introduced and coded as -4 (the code assigned 
to cDs in the canonical Trh system). 

Before illustrating in some detail the MORPHOT tool, 
it is important to explain what would be the ideal perfor- 
mance that we have tried to achieve. This limit is obviously 
represented by the intrinsic uncertainty of the morphological 
classifications provided by experienced human classifiers. 

In order to quantify this ideal target, we have first col- 
lected SDSS g-band images of 163 galaxies in the redshift 
range 0.005-0.015, with FWHM^2 arcsec and Trh classi- 
fication code given in the RC3. In this preliminar sample 
the fraction of early-type galaxies turned out to be small 
compared with that in clusters. Therefore, we decided to 
include in this sample 70 more SDSS images of early-type 
galaxies obeying the same previous criteria about FWHM 
limit and morphological availability on the RC3, but in the 
redshift range 0.015-0.03. It is worth noting that, in spite 
of the worse image quality (average seeing and photometric 
depth) of SDSS with respect to WINCS, due to the lower 
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Table 1. Revised Hubble Type (Trh) and MORPHOT Type 
(Tm) codes 



Code 


Trh 


Tm 


-O 




CJJ 


-0 


rLi 




A 

-4 


CU 


/Qn 


Q 

-o 


Qn— 


Qn— 


-2 


SO 


SO 


-1 


S0+ 


S0+ 





so/a 


so/a 


1 


Sa 


Sa 


2 


Sab 


Sab 


3 


Sb 


Sb 


4 


Sbc 


Sbc 


5 


Sc 


Sc 


6 


Scd 


Scd 


7 


Sd 


Sd 


8 


Sdm 


Sdm 


9 


Sm 


Sm 


10 


Im 


Im 


11 


cl* 


cl* 



* cE and cl are the compact elliptical and irregular galax- 
ies, respectively. 



range of redshift (z^0.03), most SDSS images of galaxies in 
our sample have visual classification accuracy (spatial res- 
olution in kpc) at least comparable to that of the WINGS 
survey (0.03<z<0.07; see Figure 8 in Paper I). 

The 233 galaxies in the final sample have been inde- 
pendently classified by two of us (AD and GF) using the 
Trh code adopted in RC3. The top panels of Figure [2] show 
the paired comparisons of the classifications from AD, GF 
and RC3, with the number of galaxies reported in each bin. 
Note that, since in the RC3 and AD databases very uncer- 
tain classification has been assigned to 42 and 8 galaxies, 
respectively (with 5 galaxies having uncertain morphology 
from both RC3 and AD), the AD-GF, GF-RC3 and AD-RC3 
comparisons just rely on 225, 191 and 188 common galaxies, 
respectively. Note also that compact galaxies, both ellipti- 
cals {Trh =-6) and irregulars {Trh = 11), are not present in 
the selected sample of RC3 galaxies. The histograms of the 
differences for each pair of classifiers are reported in the bot- 
tom panels of Figure |2] while the main statistical quantities 
of these differences are reported in the first three rows of 
Table |2] The worse performances of the RC3 clas- 
sifications with respect to those given by AD end 
GF can be explained because of the very nature 
of the RC3 data, which mainly result from com- 
pilation and statistical homogenization of different 
(mostly inhomogeneous) data sources. The fourth and 
fifth rows of Table [2] report the same quantities relative to 
the comparisons of the morphological type estimates given 
by two of us (GF and WJC) for the clusters Abell 1643 and 
Abell 1878 (z^0.20 and z^O.25, respectively; ground based, 
very good seing imaging taken at the NOT) and for the 
clusters CL 0024+16 and CL 0939+47 (z-0.39 and z-0.41, 
respectively; WFPC2 imaging from the MORPHS collabora- 
tion; Small et al.|1997 ). These comparisons are illustrated in 



Section |4j From Table[2]the visual morphological classifica- 
tions turn out not to be biased among each other, the largest 
average displacement in the table being less than AT ^0.5. 
Instead, both the rms and the fractions of absolute differ- 
ences less than 1 , 2 and 3 times Trh codes turn out to share 
relatively wide ranges {a at from 1.2 to 2.4; |AT| ^1 from 
-0.53 to -0.84; |AT| ^2 from -0.79 to -0.96; |AT| ^3 
from —0.91 to —0.99). It is interesting to note that similar 
uncertainties on the visual classifications and similar wide 
ranges in the statistical quantities of the differences were 
found by the MORPHS collaboration in their morphologi- 
cal catalog of 1857 cluster galaxies at z— 0.5, observed with 
WFPC2@HST and classified by 4 different human classifiers 
dSmail et al1|1997| their Figure 1). 



3 THE MORPHOT TOOL 

Figure [3] shows the flow-chart of MORPHOT. The top and 
bottom parts of the figure illustrate the calibration and ap- 
plication stages, respectively. Each stage must be read fol- 
lowing the direction of the big arrow on the left side. In 
particular, in the calibration stage, the visual estimates Ty 
(in the MORPHOT system Tm) are obtained for two sam- 
ples, each one including —1000 galaxies, extracted with the 
same random ctiteria from the WINGS imaging. The first 
one will be used as a calibration sample for the tool, while 
the second one will be employed in Section |4] as a test sam- 
ple in order to assess the performances of MORPHOT. For 
each galaxy in the calibration sample, (i) the global quan- 
tities: size (R), signal-to-noise ratio {S/N) and ellipticity 
{s) are recorded; (ii) 20 image-based, numerical diagnostics 
of morphology {Di, i=l,...,20) are defined and their values 
are evaluated. The calibration sample is used to gauge how 
the diagnostics Di depend on Ty and on the global quan- 
tities. This allows us to produce a semi-analytical estima- 
tor which combines the most effective diagnostics through a 



Maximum Likelihood technique (ML; see Section 3.3.1 and 
Appendix [b]). The same sample is also used as a training set 
for a Neural Network machine (NN; see Section 3.3.2 and 
Appendix [C]), in which the global quantities (R, S/N and 
e) and the diagnostics Di are the input quantities and the 
visual codes Ty (in the Tm system of Table [T]) are the tar- 
gets. Finally, the NN and the ML estimators are combined 
to produce the final MORPHOT estimator Tm- 

In the following sub-sections of the present Section the 
various steps of the MORPHOT tool are described in detail. 



3.1 The calibration sample of WINGS galaxies 

In the framework of the WINGS project, we have devised 
the multi-object, automatic surface photometry tool GAS- 
PHOT dPignatelh eTaL][2006l fP'Onofrio et al.||2011| ). This 
tool has been used to perform the detailed surface photom- 
etry of 42297 galaxies in the WINGS clusters for which 



Fasano et al. (2000 their Figure 2). Comments about the 
last row in the Table are given at the beginning of 



SExtractor (Bertin and Arnouts 1996) found more than 
200 (300) contiguous pixels (threshold area: Athr) 
brighter than 1.25 (1.07) times the rms per pixel 
of the background ((Jbkg) for those images obtained 
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Figure 2. The top panels show the paired comparisons of the classifications from AD, GF and RC3, with the number of galaxies reported 
in each bin. The bottom panels report the histograms of the differences for each pair of classifiers 



Table 2. Comparisons among visual morphological classifications 



Comparison A^^^^ z Telescope < AT > ctat |AT| ^1 |AT| ^2 |Ar| ^3 



AD-GF 


225 


^0.03 


SDSS 


0.076 


1.257 


0.836 


0.960 


0.982 


GF-RC3 


191 


^0.03 


SDSS 


-0.554 


2.374 


0.529 


0.796 


0.932 


AD-RC3 


188 


^0.03 


SDSS 


-0.425 


2.272 


0.569 


0.787 


0.910 


GF-WJC 


67 


-^0.2 


NOT 


-0.242 


1.348 


0.727 


0.909 


0.985 


GF-WJC 


207 


-0.5 


HST 


0.043 


1.479 


0.773 


0.928 


0.976 


GF-GF 


136 


0.04^0.07 


INT+MPG 


-0.072 


1.158 


0.940 


0.976 


0.994 



with the WFC@INT (WFI@ESO) Q (- /iv=25.7 for 
the WINGS survey). 

With the aim of providing a sample of galaxies suit- 
able to calibrate MORPHOT, optimizing its performances 
for WINGS, we decided to randomly extract 16 galaxies per 
cluster from the WINGS-GASPHOT catalogs, taking care 
of putting two galaxies in each one of the 8 bins of appar- 
ent V magnitude defined as follows: V^15, 15<V^16, ... 
, 20<V^21, V>21. In this way we gathered 1216 WINGS 
galaxies sampling uniformly the whole range of magnitudes 



^ The different thresholds cr^fc^ and number of contiguous pix- 
els take into account that we are using cameras with different 
pixelsize. See jVarela et al. ( 2009| for details. 



of the GASPHOT-WINGS galaxy sample (see Section [1} 
and spanning the whole range of observing conditions (back- 
ground noise and FWHM) of the WINGS optical imaging 
(see Varela et al.|f2QQ9 ). We decided to remove from this 
sample those galaxies too close to the edges of the CCDs 
and/or the very peculiar objects (on-going mergers or quite 
ill-shapen galaxies). After that, we are left with a final cali- 
bration sample of 926 galaxies. All these galaxies have been 
visually classified by GF according to the Tm code shown in 
Table [T] 

It is worth recalling here that the photometry of the 
WINGS optical survey has been performed on images in 
which large galaxies and halos of bright stars have been 
removed after modeling them with elliptical isophotes (see 
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MORPHOT flow-chart 



calibration sample 

-1000 WINGS galaxies 
(16 per cluster, randomly 
selected in 8 bins of mag) 



< 

a: 
cn 



visual type (TJ 



< 



size (Rgo) and 
(S/N) of galaxies in 
the calib. sample 



definition of diagnostics 
Di (i=1,...,20 + ellipticity£) 
(see Appendix A) 



post-stamp 
images 



image 
pre-processing 



Neural Network 

training 
(see Appendix C) 




diagnostics (Dj+e) 
of galaxies in the 
calib. sample 

L 



Max. Likelihood 

calibration 
(see Appendix B) 



Neural Network 
estimates 



Max. Likekihood 
estimates 




final MORPHOT 
estimates 



diagnostics (Dj), size (Rgg) and (S/N) 
of galaxies in the WINGS clusters 



1 



image 
pre-processing 



post-stamp 
images 



1 



WINGS cluster 
image 




Figure 3. Flow-chart of MORPHOT 



Galaxy morphology in WINGS 7 



Varela et al.|2009 ). This allowed us to perform a careful sub- 
traction of the background and estimation of its rms {(Jbkg) 
even for small galaxies embedded in the halo of the bright- 
est cluster galaxies. Therefore, the WINGS optical catalogs 
provide, for all galaxies, robust SExtractor determination of 
the ellipticity (ssex) and of the above mentioned threshold 
area Athr- 

It is also worth mentioning that each individual galaxy 
is recorded within a square, odd-sized frame of side 3xamaj, 
where amaj corresponds to the semi-major axis (in pixels) 
of the ellipse with area Athr and ellipticity ssex- 



3.1.1 Preliminar image processing 

Before running the core-tool of morphological type estima- 
tion, for each post-stamp galaxy image of a given sample, 
MORPHOT automatically performs the refinement of the 
local background subtraction and the galaxy re-centering. In 
particular, the central pixel of each galaxy image is made to 
be coincident with the intensity peak or with the distance- 
averaged intensity (bari) center, depending on whether the 
galaxy shows a well defined, dominant light peak (regular 
shape) or an irregular structure with several local peaks. 
Moreover, a preliminar processing of the post-stamp galaxy 
image is performed, which produces two ancillary (and tem- 
porary) frames Fc and Fs. 

In the frame Fc (C: Clean), the possible spurious fea- 
tures (ghosts) and/or those objects (both stars and galax- 
ies) different from the galaxy under analysis are removed by 
comparing each other the original (Fo)and the 180° rotated 
(Fiso) images. In particular, if at a given pixel position the 
difference Fq-Fiso is greater than n times the rms of the 
pixel values of Fiso over a box of side = 2xFWHM around 
the same position, the pixel value in Fq is replaced by the 
corresponding value in Fiso- We have empirically verified 
that, in our redshift range and with our instrumental set, 
using n=3 allows us to remove satisfactorily most of the un- 
wanted objects without either changing the statistical prop- 
erties of Fq (image texture), or fading the interesting galaxy 
features, like spiral arms, bars, rings and HII regions. 

The frame Fs (S: Smooth/Symmetric) is obtained in 
two steps: first the symmetrization is achieved by averaging 
(pixel by pixel) Fc with its 180° rotated version; then, the 
median (3x3) and adaptive (Richter et al.|19"9T max. block 
size=ll) filters are applied to the symmetrized frame. 

It is worth stressing that the evaluation of the mor- 
phological diagnostics (see following Sections) is performed 
on either Fc or Fs, depending on the particular diagnos- 
tic. In general, those more specifically linked to the global 
properties of galaxies (i.e.: different kinds of concentration, 
etc..) are evaluated on Fs, while those dealing with pixel- 
scale structures, local features (dumpiness, diskyness, etc..) 
and simmetry are evaluated on Fc. The cleaned image Fc 
is also used to determine the total intensity (It, in ADUs) 
and the final, global geometrical parameters (ellipticity e 
and position angle 0) of the galaxy. These are used in turn 
to produce a model image Fm from the elliptical apertures 
intensity profile of Fc, to determine the equivalent radii (in 
pixels) enclosing 80% of the total galaxy light (Rso) and 
to compute the average signal-to-noise ratio of the galaxy: 
S/N=It/ {(Jbkg Athr)- The global quantities e, Rso and S/N 



are used in the calibration procedure of diagnostics (see Sec- 
tion 



3.2) 



3.2 Morphological diagnostics 

Our approach to the automation of morphological classifi- 
cation is a fully empirical one. We do not try to identify an 
orthogonal set of few, independent and physically pregnant 
morphological indicators (diagnostics, hereafter), as in the 
case of the CAS parameter set (Concentration, Asymmetry 
and dumpiness; Consel ice|2003 ) or in the papers by van der 
Wl([2008) and Scar lata et al.|(2007|). Rather, we prefer to 



bet on a large number of diagnostics, no matter if in some 
cases they are similar among each other, since we postulate 
that each one of them could potentially be sensitive to some 
particular morphological characteristic and/or feature of the 
galaxies. In other words, we decided not to throw out any- 
thing a priori and to defer to a later stage the possiblity of 
giving up some of the diagnostics. As well, we do not try 
to select the most significant diagnostics by means of statis- 
tical techniques, like for instance the Principal Component 



Analysis (PCA, see Section 3.3). We just test each diag- 
nostic on the field (the test sample, see Section [4]), checking 
whether its addition to the previous (smaller) set of diagnos- 
tics improves the tool's performances (see Section [3.3.1 and 
Figure [s]). Lastly, our diagnostics are not necessarily defined 
(and conceived) to be independent of the image parameters 
(photometric depth, noise, pixel size, seeing, etc.). In fact, 
at least for the time being, we aim to apply MORPHOT 
just to the WINGS imaging, deferring to a later time the 
release of a more generally usable version of the tool, where 
the definition of the diagnostics will be as most as possible 
independent on the observing material. More explicitly, the 
calibration of the diagnostics we describe in the next Sec- 
tions is performed on the WINGS calibration sample defined 
in Section [3T] and holds good just for WINGS-like data. For 
now, applying MORPHOT to imaging data different from 
WINGS would actually imply a re-calibration of the diag- 
nostics on the new data set. 

Up to now we have empirically introduced (and tested) 
20 diagnostics. Some of them are not conceptually origi- 
nal, but are usually more simply defined with respect to the 
similar indicators already present in the literature. Again, 
we prefer to test a large number of rough diagnostics that 
are rapidly evaluated, rather than a small set of carefully 
calibrated (but sometimes hard to compute and not neces- 
sarily more efficient) indicators. It is also worth noting that 
our set of diagnostics is actually open, meaning that addi- 
tionally devised diagnostics (like, for instance, the spirality 
analyzers from [Naim et al.|1995 Shamir|2011 ) can be intro- 
duced without changing the structure of the tool. The only 
limitation we pose for the new diagnostics is that they have 
to be image-based, thus excluding color- and spectroscopy- 
based quantities. This is because we think that, in order 
to have an unbiased picture of the evolution of galaxies in 
clusters, the information on morphology and stellar popu- 
lation should be kept separate. In Appendix [A] we present 
in some detail the definition and the meaning of the 20 di- 
agnostics Di (i=l,...,20). Here we just mention that many 
of them turn out to be correlated, sometimes strongly so. 
This is shown in Figure^ for the calibration sample and it 
is somehow expected due to our empirical approach. 
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Figure 4. Mutual correlations among the MORPHOT diagnostics defined in Appendix [A] 



S.2.1 Diagnostic dependences 

Figure |5] illustrates how the 20 morphological diagnostics 
defined in Appendix [A] correlate with the visual morpho- 
logical type for the 926 galaxies of the calibration sample. 
For the sake of clarity the Y-labels are omitted in the figure 
and the visual morphological types (integer values) of the 
calibration galaxies in the southern clusters (304 objects) 
are shifted by 0.5 upward. 

Moreover, in Figure [5] the diagnostics Di, D7, Dg, Dq, 
Dio, D12, Dis, Di9 and D20 are plotted in logarithmic scale. 
It is evident from the figure that many diagnostics have 



quite similar (average) behaviour as a function of the vi- 
sual morphological type. Still, as explained before, we pos- 
tulate that even slight differences among diagnostics could 
in principle help to disentangle different morphological fea- 
tures and we defer the final decision about the diagnostics to 
be retained to the comparison of the results with the visual 
classifications. For the moment, it is worth emphasizing the 
importance of the diagnostics Dio and Du in disentangling 
ellipticals from SO galaxies. Figure [6] shows that the distri- 
butions of these diagnostics for the two morphological types 
are quite apart. The relative scarcity of non-disky (face-on) 
SOs in the figure is likely attributable to some cases of face-on 
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Visual Type (T,) (Djo + Dn)/2 



Figure 5. The MORPHOT diagnostics defined in Ap- 
pendix vs. visual morphological type for the calibration 
sample. 



Figure 6. Distribution of the average value of the diskiness diag- 
nostics Dio and Dn (see Appendix [a]) for the (visually clas- 
sified) elliptical and SO galaxies of the calibration sample 



SOs which have been visually misclassified as ellipticals. To 
this concern we note that, with the spatial resolution typ- 
ical of the WINGS survey, to entirely remove this kind of 
miss-classification turns out to be almost impossible, even 
for visual classifications. Actually, [Capaccioli et aL (1991) 
have shown that to distinguish face-on SOs from ellipticals 
could be a difficult task even for very well resolved galaxies 
(see Cappellari et al.|2011| for a more radical point of view). 

As already stated before, since our aim is to pro- 
vide morphological classifications of WINGS galaxies, our 
diagnostics D^ are not conceived to be independent on 
instrumental and observing parameters (pixelsize, seeing, 
S/N ^ etc.). Figure |7| illustrates how the Di depend on 
\og{RsQ/FWHM) and on log(5'/A^). 

Instead, no significant dependence of the diagnostics on 
the apparent ellipticity e has been found. To this concern, 
before describing the techniques we used to extract from our 
diagnostics univocal estimates of the morphological type, it 
is worth mentioning that, from now on, we formally include 
£ into the set of diagnostics. Therefore, their total number 
is hereafter assumed to be of 21 (1)^,1=1, ..,20+s). 



3.3 Combining the diagnostics 

Having defined the diagnostics Di and tested their depen- 
dence on both the visual morphological type {Ty) and the 
global quantities Rso and S/N ^ we are left with the diffi- 
cult task of simultaneously exploiting their capabilities, in 
order to improve as much as possible the final effectiveness 
of the tool in recognizing the morphology of galaxies. In 
other words, we must combine in some (smart) way the 21 
diagnostics to obtain a single, final morphological estimator. 
Although our empirical approach would drive us to use 



all the diagnostics (see Section 3.2), we first tried to iden- 
tify, through the canonical Principal Component Analysis 
(PGA), an orthogonal transformation converting our diag- 
nostics in a set of uncorrelated variables, smaller than the 
original set, but still preserving the wealth of morphologi- 
cal information contained therein. However, likely because 
our diagnostics are not normally distributed, this attempt 
turned out to be unsuccessful. In fact, running the PGA on 
our galaxy calibration sample, we just obtained two signif- 
icant eigenvectors, whose linear combination resulted in an 
extremely large scatter of the PGA morphological types with 
respect to the visual estimates. Actually, in Section 3.3.1 (see 
Figure [s]) the number of significant diagnostics is shown to 
be much larger than two. 

Returning to our empirical approach, we used two dif- 
ferent techniques, totally independent on each other, in or- 
der to obtain the above mentioned combination of the diag- 
nostics and the final, global morphological estimator. It is 
worth mentioning that both techniques produce morpholog- 
ical type estimates in one digit decimal numbers. 

3.3.1 Maximum Likelihood estimator 

As outlined at the beginning of Section [s] (see also the flow- 
chart in Figure Bl), the first technique exploits the Maxi- 



mum Likelihood (ML) statistics to combine the diagnostics. 
Goncisely, after having removed their obvious dependences 
on the galaxy size (relative to the FWHM) and signal to 
noise ratio (see Section 3.2.1 and Figure [7|, we use the de- 
pendence of diagnostics on the visual morphological type in 
the calibration sample (i.e. the 2D distributions in Figure |5] 
to estimate the probability that a given value of each di- 
agnostic could come from (be measured for) galaxies of all 
possible morphological types. Then, for a galaxy with un- 



10 G. Fasano et al 














iI m^ iI mm I mm I 

















D3 



























0.5 1 1.5 20 0.5 1 1.5 20 0.5 1 1.5 20 0.5 1 1.5 2 1.5 2 2.5 3 3.5 1.5 2 2.5 3 3.5 1.5 2 2.5 3 3.5 1.5 2 2.5 3 3.5 

Log(R3o/FWHM) Log(S/N) 
Figure 7. The MORPHOT diagnostics defined in Appendix |a] vs. \og{R8o / FW H M) and log(S'/A^). 



known morphology and known (measured) diagnostics, we 
compute the ML probabihty (product of the probabihties 
associated to the diagnostics) as a function of the morpho- 
logical type T and we assume that the Hrue^ morphological 
type of the galaxy is that providing the largest value of ML. 
From the function ML(T), we can also derive the confidence 
interval of each ML estimate of the morphological type. De- 
tails about the MORPHOT-ML technique can be found in 
Appendix |B] 

In order to determine how many diagnostics (and which 
ones) are necessary (and sufficient) to optimize the ML 
techique, we have first applied the above procedure to all 
galaxies in the calibration sample using the 21 Di one by 
one and recording the diagnostic which provides the lowest 
scatter (rms) of the differences between the ML and visual 
morphological types {Tml — Ty] hereafter ATml)- Then, 
we have repeated the procedure by adding the remaining 20 
diagnostics one by one to the firstly selected diagnostic, and 
we have again recorded the one which minimizes the above 
mentioned rms among the 20 couples of diagnostics. We it- 
erated this loop, each time adding one by one the remaining 
n diagnostics to the (21-n) already recorded, while the Di 
last. Figure [S] illustrates the result of this iteration, showing 
how the average value and the rms of the ATml distribu- 
tion vary as a function of the number Nd of diagnostics used 
to provide the Tml of galaxies in the calibration sample. 

The circled numbers in the figure refer to the corre- 
sponding diagnostic's numbers in Appendix ^ Figure [S] 
shows that the average ATml becomes nearly stable (at 
^0.24) after Nd ^9, while the rms of the ATml distri- 
bution decreases till Nd It is worth mentioning that, 
according to the F-test for significantly different variances, 
up to this value of Nd, the addition of new diagnostics sig- 
nificantly reduces the rms, at variance with the (above men- 
tioned) formal result of the PCA. From Figure [s] the diag- 




Figure 8. Average value (upper panel) and rms (lower panel) of 
the differences ATm l = {Tm L — Ty) as a function of the number of 
diagnostics (Nd) used for the ML estimator. The numbers in the 
open circles identify the diagnostics (see Appendix^ recursively 
added to the previus ones. The average value and the rms become 
nearly stable after Nd=9 and Nd=H, respectively 



nostics which turn out to be effective for the ML technique, 
sorted by decreasing effectiveness, are: D13, i^is, D9, i^i2, 
Dii, D5, £, Di7, Di4, Dio and D3. 

The outlined iterative procedure, aimed at identifying 



Galaxy morphology in WINGS 11 



the most effective diagnostics of the ML technique, is clearly 
an empirical one. For instance, we have chosen to stop 
the iterations when minimizing the rms (A/^d=11)? 
rather than the average value of ATml, since we con- 
sider the minimum ATml at Nd=^ to be just a sta- 
tistical fluctuation due to the finiteness of the test 
sample. Moreover, we cannot rule out (actually, we con- 
sider very likely) the possibility that different combinations 
of Di coud work better, giving lower values of rms{ATML)- 
Still, testing all possible combinations of the diagnostics was 
intractable and, we believe, unproductive. Thus, we assume 
that the morphological types Tml we got using the first 
eleven effective diagnostics are the best possible ML esti- 
mates for the calibration galaxies. Although lacking in a 
rigorous explanation, we believe the slightly worse 
rms performance of the tool for Nd >11 to be due 
(again) to the flniteness of the test sample, which 
might induce in the empirical ML procedure a sort 
of oversampling noise. 

The short-dashed lines in Figure [9] connect the median 
values of ATml in different bins of Ty (upper panel) and the 
corresponding rms values (lower panel), while the relevant 
global statistics of the ATm l distribution for the calibration 
sample are reported in the first row of Table |3] 



3.3.2 Neural Network estimator 

The second technique we use to combine the morphologi- 
cal diagnostics is based on the classical feed-forward multi- 
layer perceptron Neural Network (NN). Details about the 
MORPHOT-NN technique can be found in Appendix jC] 
Here we just mention that again the NN morphological 
type estimates are supplied with confidence intervals, while 
in this case (at variance with the ML technique) we use 
as input quantities of the NN machine the whole set of 
diagnostics (D^, i=l,...,20) plus the global quantities e, 
\og{R8o/ FWHM) and log(5'/A). The reason for this choice 
is explained in Appendix [C] 

The dot-dashed lines in Figure |9] (blue in the electronic 
version) connect the median values of ATnn in different 
bins of Tv (upper panel) and the corresponding rms val- 
ues (lower panel), while the relevant global statistics of the 
ATnn distribution for the calibration sample are reported 
in the second row of Table [S] Comparing these values with 
the corresponding ones relative to the ML technique (first 
row of the same table) and looking at Figure |9] we conclude 
that the performances of the NN estimator are significantly 
better than those of the ML estimator. 



3. 3.3 The final MORPHO T estimator 

As already pointed out at the beginning of Section [331 ^^e 
Maximum Likelihood and the Neural Network provide con- 
ceptually and technically different approaches to the prob- 
lem of combining the morphological diagnostics. Therefore, 
the MORPHOT estimators produced by the two techniques 
{Tml and Tnn) should be independent on each other and 
the rms of their difference should roughly be the square 
root of the sum of their variances with respect to Ty 
{rms[TML — Tnn] ^2.55). Actually, elementary numerical 
simulations show that the particular density distribution of 
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Figure 9. Upper panel: The median values of (T—Ty) in different 
bins of Ty are connected for the estimators Tml (short-dashed 
line), Tnn (dot-dashed line, blue in the electronic version) and 
Tm (long-dashed line, red in the electronic version) obtained run- 
ning MORPHOT on the calibration sample. The full (black) line 
illustrates the behaviour of the Tm estimators obtained running 
MORPHOT on the test sample (see Section^. Lower panel: The 
rms of (T — Ty) in different bins of Ty for the estimators Tml, 
Tnn cind Tm in the case of the calibration sample. The meaning 
of the different lines are as in the upper panel. The full (black) 
line refers to the Tm estimator in the case of the test sample. 



Ty makes the real rms lower than the above theoretical 
value, in agreement with the value we found in the calibra- 
tion sample (rms=2.05, see last row in Table [Sj. 

Once the mutual independence of the two estimators 
has been checked, the last step of the MORPHOT flow- 
chart (see Figure [3| is the evaluation of the final morpho- 
logical type estimator Tm , which is simply defined as the av- 
erage value of the two independent estimators: Tm={Tml + 
Tnn)/'2. Similarly, the lower and upper limits of the con- 
fidence interval of Tm are obtained by averaging the lower 
and upper confidence limits of Tml and Tnn, respectively. 

The long-dashed lines in Figure |9] (red in the electronic 
version) connect the median values of ATm in different 
bins of Ty (upper panel) and the corresponding rms val- 
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ues (lower panel), while the relevant global statistics of the 
ATm distribution for the calibration sample are reported in 
the third row of Table [S] 

Comparing these values with the corresponding ones 
relative to the ML and NN techniques (first two rows of 
the same table), we should conclude that the global per- 
formances of the NN estimator are even better than those 
of the final MORPHOT estimator, thus making convenient 
to adopt Tnn alone to optimize the performances of MOR- 
PHOT. Still, just because of the above mentioned mutual 
independence of the ML and NN estimators, we are inclined 
to believe that their combination could in any case provide 
more stable results, each technique possibly compensating 
the biases of the other one. Actually, from the upper panel 
of Figure |9] the biases of the Tm and Tnn estimators in dif- 
ferent bins of Tv have quite similar sizes, while the scatter 
of the Tm estimator along the Ty sequence (long-dashed, 
red line in the lower panel) turns out to be more stable than 
in the case of the Tnn estimator (dot-dashed, blue line in 
the same panel). Therefore, we decided to adopt Tm as final 
MORPHOT estimator. 

It is worth noting that all the estimators tend to be 
biased towards later and earlier morphological types for the 
early and late visual types, respectively. However, this is 
expected (and in some sense obvious) because of the one- 
sided error distribution of galaxies close to the limits of the 
available range of visual morphological types. 



4 TESTING THE PERFORMANCES OF 
MORPHOT 

As mentioned in Section [3] the test galaxy sample has been 
extracted from the V-band WINGS imaging using the same 



Test Sample 



(random) criteria described in Section 3.1 for the calibration 
sample. Again, we removed from the initial sample of 1216 
galaxies those objects too close to the edges of the CCDs 
and/or the very peculiar objects (on- going mergers or quite 
ill-shapen galaxies), being left with a final sample of 979 ob- 
jects, which have been visually classified by GF according 
to the Tm code exemplified in Table [l] It is worth noting 
that 136 (mostly bright) galaxies turned out to be in com- 
mon between the calibration and the test samples. Since the 
classifications of the two samples (both from GF) are inde- 
pendent from each other, these galaxies have been used to 
estimate the internal consistency of the visual classifications 
from GF. The main statistical indicators of this com- 
parison are reported in the last row of Table |2} They 
show that significant differences can be found even 
comparing among each other the morphological clas- 
sifications provided by the same human classifier, for 
the same galaxy sample, but at different times. 

For all the galaxies in the test sample, we have com- 
puted the global quantities (e, i^so and S/N) and the MOR- 
PHOT diagnostics Di and we have obtained the Maximum 
Likelihood (Tml), the Neural Network (Tnn) and the final 
(Tm) MORPHOT estimators of the morphological type. 

In the fourth row of Table [3] we report the relevant 
statistical quantities of the comparisons between visual and 
MORPHOT classifications for the test sample. 

From Tables |2] (rows 1, 4 and 5) and [s] (rows 3 and 
4) and from Figure |9] one can derive the following remark- 
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Figure 10. Dependence of the difference (Tm — Ty) on apparent 
magnitude (upper panel), threshold area (middle panel) and sig- 
nal to noise ratio (bottom panel) for the test sample. The big, 
full dots and the corresponding error bars illustrate the 
average values ond the rms in each bin. 



able conclusions: (i) for the calibration sample the scatter of 
the Tm estimator with respect to the visual classifications 
turns out to be quite comparable to (sometimes better than) 
the scatter reported in Table [2] among visual types provided 
by different experienced human classifiers; (ii) for the test 
galaxy sample, the above mentioned scatter just marginally 
increases with respect to the previous case, still remaining 
quite competitive with respect to the average scatter among 
visual classifications. 

Figure [To] illustrates how the differences /STm—Tm-Tv 
(average value and scatter) behave as a function of apparent 
magnitudes, threshold areas and the S/N ratios of galaxies 
in the test sample. As might be expected, the scatter of ATm 
increases (from ^^0.7 to ^2) at decreasing the threshold area 
(in pixels) and the apparent luminosity of galaxies. Instead, 
no dependence of the scatter on S/N is found. 

Figure [TT] illustrates the comparison between the vi- 
sual and the final (Tm) morphological classifications of the 
galaxies in the test sample. In this case the comparison 
is made in 'broad' bins of morphology, where the 'broad' 
classes are conventionally defined as follows: Ellipticals (E), 
for T <-4; Lenticulars (SO), for -4^ T ^0; Early-Spirals 
(SpE), for 0< T ^4; Late-Spirals (SpL), for 4< T ^7; Very 
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Table 3. Comparisons between visual morphology and MORPHOT Type for both the calibration and the test samples. The last row 
refers to the mutual comparison between Tml and Tjsin • 
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Figure 11. Comparison between visual and MORPHOT 'broad' 
morphological classes for the galaxies of the MORPHOT test 
sample. At the top of the 2D-bins the percentages of the visual 
classes (Es, SOs, early spirals [SpE], late spirals [SpL] and irregu- 
lars) falling in different bins of the MORPHOT classification are 
reported. Similarly, the percentages of the MORPHOT classes 
falling in different bins of the visual classification are reported in 
the right side of the 2D-bins. Finally, on the top (columns) and on 
the right (rows) of the plot, we report the total number of galax- 
ies in each 'broad' class of the visual and MORPHOT estimates, 
respectively. 



Late-Spirals and Irregulars (Irr), for T >7. Note that we 
have included the cD galaxies in the broad class E. This is 
because MORPHOT tends to classify as cDs (T=-6) some 
among the brightest and largest ellipticals in the test sam- 
ple. Note also that we have included in the SO 'broad' class 
both the galaxies classified E/SO (T=-4) and those classified 
so/a (T=0). 

Finally, Figure [12] illustrates in more detail the results 
already shown in the previous two figures and compares, 
for different 'broad' morphological classes, the distributions 



of apparent magnitude, threshold area, signal to noise and 
axial ratios obtained from the visual (thick lines) and MOR- 
PHOT (thin lines, red in the electronic version of the 
paper) estimates of the test galaxy sample. 

Figures [TT] and [12] show that, besides mimiching the 
statistics < AT > and a at of the comparison between hu- 
man classifiers (see Tables |2] and |3| , the automatic mor- 
phological types are able to fairly reproduce the global mor- 



phological fractions of the visual types (Figure 11), as well 
as the fractions binned according to several observed quanti- 
ties (Figure 12). In particular, Figure 11 shows that, in 



spite of the wide range of Tm corresponding to each 
bin of Tv (and viceversa), the global fractions of vi- 
sually classified E and SO galaxies in the test sample 
are almost exactly reproduced by the MORPHOT 
types. 



5 APPLYING MORPHOT TO THE WINGS 
CLUSTERS 

The bottom part of Figure [3] illustrates the flow-chart rela- 
tive to the application of MORPHOT to the WINGS clus- 
ters. 

Concisely: (i) for the galaxy catalog of a given WINGS 
cluster, post-stamp frames of the galaxies to be classified are 
extracted from the original WINGS imaging; (ii) for each 
galaxy, the global quantities (Rso, S/N and e) and the di- 
agnostics Di are evaluated and the (independent) NN and 
ML estimators of the morphological type are produced, each 
one with the proper confidence interval; (iii) the final MOR- 
PHOT estimator is obtained by averaging the NN and ML 
estimators. 

According to [Varela et al.| ( |2009| , the WINGS optical 



(B,V) imaging provides photometric and geometric descrip- 
tion of 400140 galaxies in 77 clusters (^5200 galaxies per 



cluster, on average). As already mentioned in Section 3.1 
for about one tenth of them (42297 galaxies: those with 
threshold area greater than 200/300 pixels for images 
from WFC@INT/WFI@ESO) t he surface photometry 
has been performed with GASPHOT ( [Pignatem et al.|2006 ). 
We use the GASPHOT- WINGS catalogs as input galaxy 
sample to perform the morphological analysis of WINGS 
galaxies with MORPHOT. 



5.1 The WINGS morphological catalogs 

The total number of galaxies in the GASPHOT catalogs is 
42,297. We removed from the sample those galaxies for which 
the Sersic index provided by GASPHOT coincides with the 
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Figure 12. The distributions of apparent magnitudes (V-band), threshold areas (pixels), signal to noise ratios and axial ratios for 
different 'broad' morphological classes (Es, SOs, early spirals [SpE], late spirals [SpL] and irregulars) in the WINGS test galaxy sample. 
The thick line histograms refer to the visual morphological types {Ty)-, while the thin line histograms (red in the electronic version 
of the paper) illustrate the distributions for the MORPHOT estimates {Tm) 



boundaries of the allowed range (0.5—8), which usually indi- 



cates that the fitting procedure was unsuccessful (Pignatelli 
et al.|2006| ). In this way we are left with 39651 galaxies in the 
fields of the WINGS clusters. For 527 galaxies (1.3% of this 
sample) MORPHOT produced unreliable results since it was 
not able to compute some of the diagnostics (fuzzy /faint ob- 
jects). The remaining 39124 galaxies have been processed by 
MORPHOT, which provides Tml, Tnn and Tm estimates 
of the morphological types, together with the corresponding 
confidence intervals T^^^ and T^^"^. 

We have very quickly checked on the WINGS imaging 
the MORPHOT classification of the bright galaxies (mostly 
brighter than and in 426 cases (^^1% of the total sam- 



ple; ^5% of the checked sample) we have modified the final 
classification since it was clearly wrong. Profiting from this 
visual inspection procedure, we have also manually added to 
the catalogs the morphological types of 799 bright galaxies 
(again mostly brighter than 1/=18), close to the borders of 
the frames or close to bright stars, which had been dis- 
carded 'a priori' by the GASPHOT and MORPHOT 
tools. After this manual intervention, the total number of 
WINGS galaxies for which we provide the final morpholog- 
ical type estimate (Tf) is 39923. Among them, 2963 are 
visual estimates (Ty). This latter number includes the 926 
galaxies of the Calibration Sample, the 979 galaxies of the 
Test Sample (136 of them turned out to be in common 
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morphological types and the visual estimates from the liter- 
ature looks satisfactory. In particular, the average value and 
the scatter of the difference {TM-TNair+) turn out to be -0.4 
and 1.4, respectively, to be compared with the corresponding 
values given in Tables [2] and [3] 
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Figure 13. Comparison between MORPHOT and literature mor- 
phological types (see the comments in the text). Left panel: num- 
ber of galaxies in the different bins of 'broad' morphological class 
for t he sample of 18 WIN GS-MORPHOT galaxies in common 
with|Fukugita et al.| ( [2007| . Right panel: the MORPHOT vs. |Nair| 
|et al.| ( |2010| ) morphological types for the 79 common galaxies. 



with the calibration sample; see Section|4|, the 426 galaxies 
whose classification has been modified after visual check (31 
of them turned out to be in common with the test sam- 
ple) and the 799 galaxies manually added to the catalogs. 
The full MORPHOT catalog of WINGS galaxies is available 
from the "Centre de Donnees Astronomiques de Strasbourg" 
(CDS) using the ViZiER Catalogue Service. Table [I] show 
the first few records of the catalog. 



5.2 External comparisons 

In order to provide an external check of the goodness of our 
automated classifications, we have searched the literature for 
visually classified galaxy samples having objects in common 
with our WINGS-MORPHOT sample. We only found three 
possible data samples that could be usable for our purpose, 
all of them concerning the S PSS galaxies: [Fukugita et al.| 
( [2007| , |Nair et al.| ( |2010| and |Lintott et al.| (|2011, Galaxy 
Zoo). By cross-matching these samples with our catalog, we 
found that the objects in common are 18, 79 and 2110, re- 
spectively. However, the potentially most sizeable compari- 
son sample (Galaxy Zoo) turns out to be practically useless, 
since it just provides the binary information Elliptical/ Spiral 
(no SO classification). Moreover, while the morphological res- 
olution of the classification system adopted by Nai r et aT] 
( 20101 is comparable (not equal) to the MORPHOT reso- 
lution, the F ukugita et al.| (|2007|) system only enables us to 
compare the 'broad' morphological classes. 

Figure [131 illustrates the comparison of the MORPHOT 



morphological types with the Fukugita et al. (2007) and 



Nair et ah] ( 12010 *) classifications (left and right panel, re- 



spectively). In both cases the MORPHOT results tend to 
be slightly shifted towards earlier types with respect to the 
visual estimates. However, in spite of the small number of 
cross-matched galaxies and of the different classification sys- 
tems, the agreement between the MORPHOT automated 



5.3 Morphological properties of the WINGS 
clusters 

This section outlines the main statistical properties of 
the galaxy morphology in the WINGS clusters. More de- 
tailed and exhaustive analyses about this topic will be pre- 
sented in a few forthcoming papers. Here we just illustrate 
some general morphological trend emerging from WINGS- 
MORPHOT catalogs. 

Figure [141 illustrates the distribution of the MORPHOT 
types (Tm) in the fields of the WINGS clusters. The con- 
tribution of the general field to the different morphological 
types has been estimated in two independent ways. First, 
we have counted galaxies in the general field of the Padova 
Millennium Galaxy and Group Catalogue (PM2GC; iCalvi 



et al.|2011 ), for which we have obtained MORPHOT classi- 
fications. Second, we have estimated the fraction of cluster 
members in the WINGS fields using the spectroscopic com- 
pleteness and membership functions derived for the WINGS 
survey b y j Cava et al.| (|2009). In the first case (upper panel 
of Figure 14 ) we can confidently assume the PM2GC survey 
to be nearly complete down to Vum^'^S, while the selec- 
tion of the WINGS spectroscopic sample has been extended 
down to V ^20 (central panel of the figure). In this case, 
the completeness function is mainly determined by time al- 
location and fiber crowding problems. Moreover, the spec- 
troscopic WINGS survey only includes a subsample of the 
original WINGS cluster sample (48 over 77). In spite of 
these differences, the distributions of the MORPHOT types 
in the WINGS clusters, obtained applying to the WINGS- 
MORPHOT catalogs the two (independent) statistical cor- 
rections for membership, turn out to be remarkably consis- 
tent (see the bottom panel in the same figure). The consis- 
tency is confirmed even if we use Vum^^S also for the mem- 
bership correction based on the spectroscopic completeness. 

Adopting the 'broad' morphological classes convention- 
ally defined in Section ^ we find that ellipticals, SOs and 
spiral galaxies constitute ^33%, 44% and 23% of the whole 
galaxy population in the WINGS clusters. It is worth to 
note that these morphological fractions are slightly different 
from those found in |Poggianti et al." (2009). The discrepancy 
mostly concerns the E/SO ratio and is due to the combined 
effects of two factors: (i) the limit of 0.6xi?2oo adopted for 



the clustercentric distance in Poggianti et al. (2009); (ii) the 
behaviour of the E/SO ratio as a function of the cluster- 
centric distance. Figure illustrates the point, showing the 
morphological fractions as a function of the clustercentric 
distance before and after correction for field contamination 
(right and left panels, respectively). Note the prevalence of 
Es in the inner cluster regions, which is responsible for the 
different E/SO ratio found in Poggia nti et al.| ( [200 9^. Note 
also that the field correction does not infiuence this ratio 
up to '-^0.5xi?2oo, mostly operating on the SO/Early-Spirals 
fraction in the external part of the clusters. 

Figure ^] shows the distributions of the projected el- 
lipticities in the WINGS-MORPHOT catalogs for different 
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Table 4. Sample rows of the MORPHOT WINGS catalog 



WINGSJD Cluster Tml TV^^ T^J- T^n T^!.- TV}-- Tm TIJ- T^-- Ty Tp 



WINGS J103833. 76-085623. 3 


A1069 


3.1 


-0.3 


4.6 


5.9 


1.6 


8.6 


4.5 


2.7 


6.3 




4.5 


WINGSJ103834. 09-085719. 2 


A1069 


4.3 


3.0 


5.2 


4.1 


0.5 


7.7 


4.2 


3.7 


4.7 




4.2 


WINGS J103834. 13-085030.4 


A1069 


-5.0 


-5.0 


-4.7 


-5.1 


-6.0 


-4.4 


-5.0 


-5.5 


-4.5 




-5.0 


WINGSJ103835. 85-084941.0 


A1069 


1.8 


-0.1 


3.7 


-0.2 


-3.5 


3.9 


0.8 


-0.5 


2.1 




0.8 


WINGS J103835. 89-085031. 5 


A1069 


-1.5 


-3.5 


0.5 


-3.2 


-5.1 


-1.4 


-2.3 


-3.4 


-1.2 




-2.3 


WINGSJ103836. 38-083614. 8 


A1069 


3.8 


1.6 


4.2 


5.0 


0.2 


8.9 


4.4 


3.6 


5.2 




4.4 


WINGSJ103836. 95-085300. 8 


A1069 


-5.0 


-5.0 


-4.5 


-5.9 


-6.0 


-5.1 


-5.5 


-5.5 


-4.5 


-5.0 


-5.0 


WINGSJ103837.15-085753.4 


A1069 


3.2 


0.2 


5.3 


4.4 


1.1 


7.9 


3.8 


2.5 


3.5 


3.0 


3.0 


WINGS J103837.48-083717.2 


A1069 


-1.9 


-3.8 


-0.0 


-3.4 


-4.8 


-1.6 


-2.6 


-3.6 


-1.6 




-2.6 


WINGSJ103837.93-084940.5 


A1069 


-5.0 


-5.0 


-4.5 


-0.9 


-5.2 


3.7 


-3.0 


-5.0 


-1.0 




-3.0 


WINGS J103838. 65-084938. 3 


A1069 


-2.4 


-3.5 


-1.1 


0.5 


-3.2 


4.8 


-0.9 


-2.4 


0.6 




-0.9 


WINGSJ103839. 63-084742. 6 


A1069 


2.0 


-2.3 


2.7 


-1.4 


-3.9 


0.9 


0.3 


-1.5 


2.1 




0.3 


WINGS J103840. 54-085041.1 


A1069 


-3.0 


-4.6 


0.8 


-4.0 


-5.9 


-1.4 


-3.5 


-4.3 


-2.7 




-3.5 


WINGSJ103840. 85-085046. 5 


A1069 


-4.9 


-3.8 


-2.3 


-5.3 


-6.0 


-4.5 


-5.1 


-5.6 


-4.6 




-5.1 


WINGSJ103841.43-085521.6 


A1069 


-3.7 


1.5 


3.3 


6.0 


0.9 


9.4 


1.1 


-3.8 


6.0 




1.1 


WINGSJ103841.63-085528.6 


A1069 


-4.8 


-3.5 


-2.7 


-3.5 


-5.5 


-1.4 


-4.2 


-4.9 


-3.5 




-4.2 


WINGSJ103842.12-083557.8 


A1069 


-5.0 


-5.0 


-4.1 


0.5 


-4.5 


6.0 


-2.3 


-5.0 


0.4 




-2.3 


WINGSJ103843.03-085602.8 


A1069 


-5.0 


-5.0 


-4.3 


-4.9 


-5.8 


-3.8 


-5.0 


-5.5 


-4.5 




-5.0 


WINGSJ103844.19-085609.3 


A1069 


0.5 


-2.1 


0.7 


-1.9 


-3.8 


0.4 


-0.7 


-2.0 


0.6 




-0.7 


WINGSJ103844.58-084601.1 


A1069 


4.5 


3.4 


5.4 


5.3 


1.2 


8.4 


4.9 


4.3 


5.5 




4.9 


WINGSJ103845.53-085341.2 


A1069 


3.1 


1.4 


4.0 


-3.1 


-5.3 


-0.5 


0.0 


-3.1 


3.1 




0.0 


WINGSJ103845. 73-084021. 9 


A1069 


0.3 


-1.6 


0.6 


-3.8 


-6.0 


1.9 


-1.7 


-3.8 


0.4 




-1.7 


WINGSJ103846.39-084530.6 


A1069 


-3.9 


-4.8 


-0.9 


-2.3 


-4.8 


0.7 


-3.1 


-4.5 


-3.5 


-4.0 


-4.0 


WINGSJ103846.47-084226.0 


A1069 


-0.9 


-2.6 


0.8 


0.7 


-1.7 


2.9 


-0.1 


-1.2 


1.0 




-0.1 


WINGSJ103847.59-084631.3 


A1069 


-1.4 


-2.1 


-0.2 


-5.5 


-6.0 


-4.3 


-3.5 


-5.6 


-1.4 




-3.5 


WINGSJ103848. 08-085044.8 


A1069 


-5.0 


-5.0 


-4.5 


-3.9 


-4.6 


-2.9 


-4.4 


-5.0 


-3.8 




-4.4 


WINGSJ103848.30-084259.9 


A1069 


-0.1 


-2.2 


1.1 


0.3 


-1.6 


2.3 


0.1 


-0.5 


0.7 




0.1 


WINGSJ103850.17-085336.6 


A1069 


-2.7 


-0.5 


2.2 


2.3 


-1.5 


6.3 


-0.2 


-2.7 


2.3 




-0.2 


WINGSJ103850.35-084804.5 


A1069 


-2.8 


-3.3 


-0.8 


-0.5 


-4.0 


3.8 


-1.7 


-3.0 


-0.4 




-1.7 



'broad' morphological classes. In this case we do not try to 
correct for cluster membership, since we just aim to check 
the plausibility of the distributions, also in comparison with 
the literature. To this concern, the ellipticity distribution 
of elliptical galaxies in Figure p^G] tu rns out to be in perfect 
agreement with Fasano et al. ( |1991[ see also Vulcani et al. 



2011), while for SOs and Early-Spirals the peaks of the dis- 



tributions are slightly shifted towards lower values of the 
ellipticities with respect to the corresponding distributions 
in Fasano et al. ( 1993| . However, this is not surprising, since 
our flattenings come from global, single component (Sersic) 
fitting of the galaxy image (GASPHOT), while the axial 
ratios in Fasano et al.| ( |1993| refer to the outer isophotes 
(essentially the disk components). The peculiar ellipticity 
distribution of Late-Spirals is due to the inclusion in this 
'broad' class of the irregular objects, which could be intrin- 
sically less flattened than disk galaxies. 

show the distributions of the (B-V) color for 



Figure 
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the spectroscopically confirmed members in the WINGS- 
MORPHOT catalogs and for the different 'broad' morpho- 
logical classes. Note the remarkable similarity and the small, 
but statistically significant shift between the distributions of 
E and SO galaxies. Note also the bi-modal color distribution 
of the Late-Spirals(+Irr) galaxies. 

Finally, in Figure we present the distribution of the 
Sersic index n for the different 'broad' morphological classes. 
Again, in this case we do not try to correct for cluster 
membership. The Sersic indices come from our WINGS- 
GASPHOT catalogs. As already mentioned in Section [l] 



even though a correlation between these n and the morpho- 
logical type exists, it is weak and it shows a high degree of 
degenaracy, especially for early-type galaxies (see the distri- 
butions for bright and faint Es). 



6 SUMMARY 

In this paper we have presented the morphological classifi- 
cation of ^40000 galaxies in the fields of 76 nearby clusters 
from the WINGS optical (V-band) survey. The morphologi- 
cal types have been estimated automatically, using the pur- 
posely devised tool MORPHOT, whose description takes up 
a substantial part of the paper. It combines a large set (21) 
of diagnostics, easily computable from the digital cutouts of 
galaxies, producing two different estimates of the morpho- 
logical type based on: (i) a semi-analytical Maximum Like- 
lihood technique; (ii) a Neural Network machine. The final, 
averaged estimator has been tested over a sample of ^1000 
visually classified WINGS galaxies, proving to be almost as 
effective as the 'eyeball' estimates themselves. In particular, 
at variance with most existing tools for automatic morpho- 
logical classification of galaxies, MORPHOT has been shown 
to be able to distinguish between ellipticals and SO galaxies 
with unprecedented accuracy. Even though its basic method- 
ology is robust for any set of digital images of similar spatial 
resolution and dynamic range, MORPHOT is presently cal- 
ibrated and fine-tuned to provide reliable morphologies of 
WINGS galaxies alone. Adjustments of the calibration are 
required (and are actually in progress) to make the tool more 
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Figure 14. Distribution of MORPHOT types in the WINGS 
clusters. Upper panel: The distributions for the uncorrected 
WINGS sam ple, for the general field sample (from PM2GC;|CaM] 
|et aI.|20Tl| and for the field-corrected WINGS sample up to the 
apparent magnitude V=18 are illustrated by dashed, dot-dashed 
and full line histograms, respectively. Middle panel: The distribu- 
tions for the uncorrected WINGS sample, for the spectroscopic 
WINGS members and for the WINGS sample corrected with the 
spectroscopic copleteness and membership functions derived from 
[Cava et al .| (|2009| are illustrated by dashed, dot-dashed and full 
line histograms, respectively. In this case a deeper magnitude limit 
(V=20) has been used (see text). Bottom panel: Comparison be- 
tween the field-corrected and the spectroscopy-corrected distribu- 
tion of the MORPHOT types in the WINGS clusters. In spite of 
the different magnitude limits adopted for the statistical correc- 
tions, the two distributions turn out to be quite consistent each 
one another. The same happens if we adopt for the spectroscopic 
membership correction the same magnitude limit we used for the 
statistical filed (PM2GC) correction (y^^^=18). 
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Figure 15. Morphological fractions as a function of the clus- 
tercentric distance (in i^200 units) for the WINGS galaxies be- 
fore and after correction for field contamination. The general field 
morphological fractions in different bins of apparent V-band mag- 
nitude have been derived from the PM2GC sample ( |CaIvi et al.| 
[20111 ) using the magnitude limit V=18. 
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Figure 16. Ellipticity distributions of WINGS galaxies for dif- 
ferent 'broad' morphological classes. The Irregular galaxies have 
been included in the Late-Spirals (see comments in the text). No 
correction for field contamination has been applied. 



generally usable. The WINGS-MORPHOT catalog has been 
exploited here just to illustrate the distributions of some 
relevant photometric and structural properties of galaxies 
in the WINGS clusters. In a few forthcoming papers of the 
WINGS series, we plan to perform more detailed statisti- 
cal analyses involving the morphology of cluster galaxies. 
In particular, besides the classical morphology-density and 
morphology-clustercentric distance relations, we will exploit 



the WINGS spectroscopic information (Cava et al. 2009 



Fritz et al. 2007 2011 Hansson et al. 2011) to study how 



galaxy morphology correlates with star formation rate and 
history at different clustercentric distances. 
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APPENDIX A: THE CURRENT SET OF 
DIAGNOSTICS: DEFINITIONS 

Here we present in some detail the definition and the mean- 
ing of the 20 diagnostics Di (i=l,...,20) we devised up to 
now. The first nine of the following diagnostics are actu- 
ally already present in the literature, although sometimes 
in slightly different forms. We will refer to the original pa- 
pers for details about their definitions. Instead, the remain- 
ing eleven diagnostics are presented here for the first time. 
Hereafter, in the definition of diagnostics, we use just the 
pixels above the threshold value {2abkg) and far from the 
galaxy center more than the image FWHM. 



Al Diagnostics already present in the literature 

Di: Sersic index of the luminosity profile. 

Given the FWHM, this diagnostic is evaluated on the image 
Fs according to the prescriptions given in |Trujillo et al.| 
Section 4 therein), making use of the previously 



The diagnostic Dg is defined according to Conselice (2003 



(2001 



extracted elliptical aperture intensity profile of the galaxy 



(see Section 3.1.1 ) 



D2: Luminosity-ranked Concentration index. 

Again from the image Fs and from the elliptical aperture 
intensity profile, this diagnostic is evaluated as the fraction 
of the total intensity coming from the 30% brightest pixels; 

D3: Distance-ranked Concentration index. 

Similar to the previous one, but defined as the fraction of 
the total intensity coming from the 30% pixels closest to 
the galaxy center (in units of elliptical distances). Note that 
more elaborated versions of the Concentration indices can 
be found in [Graham et"aL] ( |2001| ), |Conselice| ( |2003| ) and[Ya^ 
mauchi et al. ( 2005 ) ; 



D4: Luminosity- ranked Gini Coefficient. 

Following Abraham et al. (2003 ) and Lotz et al." (2004) , in- 
side the square whose sides coordinates are the fraction of 
galaxy pixels and the fraction of the total counts (square 
area = 1), we define this diagnostic as the area between 
the diagonal of the square and the galaxian Lorentz curve 
(i.e.: the rank-ordered cumulative distribution function of 
the pixel counts). For this diagnostic and for the following 
one we use the image Fs; 

D5: Distance-ranked Gini Coefficient. 

Similar to the previous one, but in this case the pixels are 
ranked in ascending order of elliptical distance from the 
galaxy center; 

Dq: Second-order moment of light. 

Following Lotz et al. (2004), we define this diagnostic as 



the second-order moment of the brightest 20% pixels of the 
image Fs, normalized to the same moment computed over 
the whole galaxy area; 

D7: Asymmetry. 

For this diagnostic we use the image Fc and we adopt 
the definition given in Conselice ( 1997 ; normalized square 



counts of the difference between the original and the 180° ro- 



tated image) , with the improvements suggested in Conselice 
|et al.| ([2000, ; see their Section 3.4) about the preliminar 
image processing (careful centering) and the handling of the 
uncorrelated noise; 

Dg and Dq: dumpiness. 



see equation 2 therein), including in the sum just the pix- 
els with counts above the threshold {2abkg) and far from 
the galaxy center more than the image FWHM. The sec- 
ond dumpiness diagnostic (Dg) is similar to the previous 
one. However, in order to further enhance high-frequency 
features, in this case the model image Fm defined in Sec- 
tion |3]lT] is subtracted from the galaxy image Fc, instead 
of the gauss-smoothed version of the image itself; 



A2 New Diagnostics 

The following morphological diagnostics are presented here 
for the first time. They all are computed on the frame Fc, 
again using just the pixels above the threshold {2ahkg) and 
far from the galaxy center more than the image FWHM. 
We recall (see Section 3.1) that Fc{i,j) is a square matrix. 



whose size (N) must be an odd number. 
Dio and Dn: Diskyness. 

The difference between the galaxy image Fc and the model 
image Fm (residual image; hereafter Fr) is used to de- 
vise two new diagnostics related to the shape of the galaxy 



isophotes (a4 > =disky, 04 < =boxy; Bender and Moel- 
lenhoff||1987 ). The first diagnostic is defined as: 



Dio = {<Fm> - < Fr2 >)/ < \Fr\ > 

where < Fri > and < Fr2 > are the average counts of 
Fr in the equal-area sectors of the above defined model- 
ellipse marked respectively with 1 and 2 in the bottom panel 
of Figure 



Al 



In the formula, the quantity < \Fr\ > is a 
normalization factor representing the average value of the 
(absolute) counts of Fr over the whole model-ellipse (apart 
from the inner circle of radius=FWHM). Clearly, Dio tend 
to be positive in galaxies with disk-shaped isophotes, since 
in this case the residuals in the two sectors marked 1 tend 
to be greater (darker in the figure) than those in sectors 
marked 2. 

The second diagnostic (Dn) is defined as the correlation 
coefficient (C.C.) between the counts of the pixels within 
whatever half-part of the model-ellipse (for instance the 
bottom-half: Ypix ^Ycen) and those of the corresponding 
pixels symmetric with respect to the galaxy center. In Fig- 
ure |A2| these two quantities are plotted against each other for 
the same galaxies used in Figure |A1| i.e. a SO galaxy with 
disky isophotes (left panel) and an elliptical galaxy (right 
panel). It is evident that high values of C.C. correspond to 
strongly disk-shaped objects, while for regular (non disky) 
ellipticals the C.C. values are close to zero; 

D12: Bandwidth of Power Spectrum. 

Roughly speaking, late-type galaxies are dominated by 
structural features (spiral arms, clumps, tails, blobs, etc..) 
whose size is (much) lower than galaxy size, while ellipti- 
cal and (in general) early-type galaxies are typically domi- 
nated by a single, regular structure, whose size is comparable 
with that of the galaxies themselves. Having this in mind, 
we have devised a new morphological diagnostic defined as 
the ratio between some characteristic inverse frequency of 
the 2D power spectrum of the galaxy image (i.e. the typi- 
cal size of features) and the equivalent threshold radius of 
the galaxy (y^ Athr/71). Operatively, we estimate the char- 
acteristic size of galaxy features by processing the Fc image 
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Figure Al. Illustration of the diagnostic Dio (see text for a 
detailed explanation). 
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Figure A2. Illustration of the diagnostic Dn (see text for a 
detailed explanation). 



with the powerspec tool (option: center=yes) included in 
the IRAF stsdas-Fourier package (FFT), and computing 
on the powerspec image the equivalent radius of the area 
where the power exceeds half of the maximum FFT power. 

The next four morphological diagnostics concern the 
statistical behaviour of very local pixel properties (image 
texture; a similar approach can be found in [Moore et al.| 



2006) over the whole galaxy body. In particular, they con- 
sider the texture unities {TUs hereafter) provided by all 
the 3x3 pixel squares centered on each pixel of the frame 
Fc{i^j) for i — j =2,...,(N-1). In order to make easier the 
formalism related to these diagnostics, it is convinient to in- 
troduce the following definitions relative to each TUij (see 



Figure A3 ) 



fij = Fc{iJ); 




Figure A3. Scheme of a Texture Unit (TU^j), illustrating the 
meaning of its upper indices (see the formulation in the text). 



ft~ = Fc{i 1); f:+ = Fc{i + 1, j + 1 



Fc{i,j-iy,fl+ 



ft- = Fc{i + i,j-iy,ft, 



JlJ 
COS aiji 
cos aij3 



= Fc{i-hj + l); 

\i-ic\/dij] cosaij2 = \i + j - ic - jc\/dijV2; 

y-ijA = \i - j\/dij\/2, 



where ic = jc =(N/2+0.5) are the coordinates of the center, 
dij — Y^(i — ic^ + {j — jcY is the distance from the center 
of the pixel (i,j) and aijk (/c=l,..,4) is the angle between 
the direction (k—) ^ in the TUij of Figure A3 and the 
line connecting (i, j) with the galaxy center. 

D13: Average Concaveness. 

Again roughly speaking, the earlier the morphological type, 
the lower the unevennes of the intensity surface of the galaxy. 
Moreover, while in early-type galaxies the intensity gradient 
increases regularly toward the center over almost the whole 
galaxy body, for late-type objects, due to the presence of rel- 
atively small structural features (spiral arms, clumps, tails, 
blobs, etc.), such regular behaviour is limited to the very in- 
ner part of the galaxy (bulge). That being stated and given 
the above definitions, this new diagnostic is expressed by the 
formula: 

A^-l 4 

^ X sgn{"^cosa,juUtj~ + /f/ - 2/^,)} 



- 0.5+ 



2 X ^ 



which in fact provides the fraction of the total galaxy lu- 
minosity coming from pixels (i^j) for which is positive the 
local concaveness, computed in the corresponding TUij and 
weight- averaged according to the direction of the galaxy cen- 
ter (cos aijk). In the above formula, sgn is the sign function: 
s^n(x)=-l,0,l for x <0, x=0 and x >0, respectively. 

Di/^i Monotonicity. 

Likewise the previous diagnostic, this one too deals with 
some geometrical rule to which the intensities inside TUs 
should obey. In particular, in this case we consider the frac- 
tion of the total galaxy luminosity coming from pixels (i, j) 
for which Fc{i,j) has a monotonic behaviour in all the 
four directions of TUij illustrated in Figure \A3\ Again, the 
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greater the amount of structural features (late-type galax- 
ies), the lower the expected fraction. In formula: 



where: 

pk±_ 



^ f ij \'f^ij ^ ij 



fc+l 



sf^Mft^ - I and <5/f -=!//;- - 



J ij 



Di^i Alignment. 

In defining the next two diagnostics it is convenient to con- 
vert the pixel coordinates in the reference system of the cir- 
cularized ellipse, whose ellipticity and position angle (e, 0) 
are those previously determined (Section 3.1.1 ) and assumed 
to be the global geometrical parameters of the galaxy: 

- jc)sin6'] X ^/Y 
-ic)sin6']/7r^ 



[{3- 



ic) cos 6 + {j - 
- jc) cos6 -\- (i - 



Using this system of coordinates, the new diagnostic (^15) 
tries to quantify to which degree the local maximum inten- 
sity gradient is aligned with the galaxy center over the whole 
galaxy body. Again, the presence of small scale structures in 
late-type galaxies should imply a lower degree of alignment 
of the local maximum gradient toward the galaxy center 
with respect to early-type galaxies. To quantify such degree 
of alignment, the points in the space {x,y, z) corresponding 
to the nine pixels of each TUij have been linearly interpo- 
lated (x^) by the function: Zij{x, y) — aij + hijX+djy. Then, 
for each TUij , we have computed the angle between two 
planes, both passing through the point [xij , yij , Zij {xij , yij)] 
and parallel to the z axis. The first plane contains the galaxy 
center (origin of the new coordinates system) , while the sec- 
ond one is parallel to the line of maximum intensity gradient 
in the previous linear interpolation. Finally, the cosines of 
the angles have been weight- averaged over the whole 
frame, according to the local intensity. In formula: 



, Zij X COS0ij 



2 



where the cosines: 



turn out to be positive if the function Zij(x,y) increases 
toward the galaxy center at the point (i,j), while they are 
negative in the opposite case. 

Die: Intercept Angle. 

In defining the previous diagnostic we have introduced the 
local planes Zij{x,y) interpolating the nine points {x,y,z) 
of each TUij. We noted that in early-type galaxies the line 
of maximum gradient of these planes should be (on aver- 
age) more oriented towards the center than in the case of 
late-type galaxies and that such alignment should be more 
and more pronounced at increasing the intensity and de- 
creasing the distance from the galaxy center. Here we note 
that, in addition, the average values of the intercept (aij) 
between these planes and the z axis should be higher in 
early-type than in late-type galaxies, where the irregular 



intensity surface should make almost randomly distributed 
both the orientations and the intercept levels of the planes, 
especially in the intermediate and outer regions of galaxies. 
In order to quantify such average intercept level avoiding 
computational divergences, we actually prefer to deal with 
the cosines of the angles {/3ij) between the z axis and the 
stright lines connecting the points (xij, yij, 0) and (a^j, 0, 0). 
Since the values of cos /3ij obviously depend on the units used 
to measure intensities and radii, it is necessary to normal- 
ize both quantities. We decided to normalize the intensities 



(intensity of the galaxy center) and the 



at Zc - Fc(ic 

radii at: Rthr — \/ Athr/T^ (equivalent radius of the thresh- 
old area). After some algebra, the new diagnostic, that is 
the normalized, average cos {3 is defined as follows: 



where: dij — y^^^IT^*^"^ circularized distance from the 

galaxy center, K = Zc/Rthr is the normalization factor and 
Wij — Zij j dij are the weighting factors, which increase at 
increasing the local intensity and at decreasing the distance 
from the galaxy center. 

Di7 — D20'. Intensity distribution moments. 

The last four diagnostics simply concern the moments of the 
intensity distribution within the galaxy frame. In particular, 
we consider the median {D17), the standard deviation (Dis), 
the skewness {D19) and the kurtosis {D20) of the distribution 
of the pixel intensities, normalized to their average value. 



APPENDIX B: THE MAXIMUM LIKELIHOOD 
TECHNIQUE 



The flow-chart in Figure [BT] should help the reader in this 
section, leading her/him through the various steps of this 
semi- analytical technique. As already mentioned in Sec- 
tion [3]2T] besided obviously depending on the visual mor- 
phological type (Figure|5|, the diagnostics Di turn out to de- 
pend on both the seeing normalized size [log^Rgo / FW H M)] 
and the signal-to-noise ratio [log(5'/A^)] (Figure |7|. There- 
fore, prior to use them to gauge the morphology of galaxies, 
we have to remove these two dependences, thus picking out 
the net dependence of each diagnostic on the morphological 
type. 

To this aim the whole range of morphological types in 
the calibration sample has been divided into eight, nearly 
omogeneous intervals (see the plots on the right side of Fig- 
ure B2). For each interval ATj (j=l,...,8), the dependences 
of all the diagnostics Di (i=l,...,21) on Rgo and S/N, have 
been represented (through least square fitting) by the fol- 
lowing second order polynomial functions: 



Dij{Rso,S/N) 



m+n=0,l,2 



' \og{R8o /FWHM)'' \og(S/Ny 



thus producing, for each Di and for each ATj, the coeffi- 
cients A^-'^ and the rms scatters Si j of the residuals with 
respect to the fits. Then, for each Di, these quantities have 
been interpolated, as a function of the morphological type, 
with Natural Cubic Sphnes [NCS)\ 
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Maximum Likelihood flow-chart 



visual type (r„), size (Rjo) 
and [S/N) of galaxies in tlie 
calibration sample 



diagnostics (Dj. r..._2o+£) 
of galaxies in the 
calibration sample 



< 

CD 



Q. 
< 



i 



for intervals of 7^ (Ty e Ar. , j = i,..,8), 
2"^^ order polynomial fit: 

D.j(R,„S/N) = 2^:;;{\ogRJ"'[\og{S/N)Y 

m+n=0,l2 

s.j=rms[D,-D.j{fit)\ 



ar(T)=Ncsr[Aru)\ 

a,(T)=NCS,[sXj)] 
(i = 1,...,21; ; = 1,...,8; m + n= 0,1,2) 



for each galaxy in the calibration sample, find the Max. 
Likelihood Type (T^l) maximizing the product of probs: 



log[n,^/(^)]' PXT) = SNF,[ADXT)] 



AD.(r) = 



try different 
sub-sets of Dj 
adopting the 

one that 
minimizes the 

quantity: 
rms (Ty-T^^) 
in the 
calibration 
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diagnostics (DJ, 
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of galaxies in the 
WINGS clusters 



Maximum 
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estimates 



Figure Bl. Flow-chart of the MORPHOT Maximum LikeUhood tool 
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Figure B3. Distribution of the deviates of the morphological 
diagnostics of the calibration galaxy sample with respect to the 
polynomial+spline fitting (details in the text) performed in order 
to remove their dependence on the size (Rgo) and signal to noise 
ratio (S/N). 



An example of this procedure is shown in the plots on the 
left side of Figure |B2| which is referred to the diagnostic 
Di3, defined in Appendix [A] Because of the evident dis- 
continuity between cD and E galaxies, the first bin (ATi) 
has been excluded from the spline fitting. Such discontinu- 
ity is actually found also for the coefficients A^''^ of the 
other diagnostics and reflects their peculiar behaviour for 
Tm——Q (see Figure [5]). It is likely telling us once again that 
cD galaxies are quite apart from Es and, in general, from all 
'normal' galaxies, thus suggesting they are the product of a 
peculiar evolutionary path (|Fasano et al. 2010). The right 
side of Figure B2 shows (again in the case of D13) the resid- 
uals of the above mentioned polynomial fits with respect 
to both \og{RsQ / FW H M) (upper panels inside each /STj 
heavy box) and log(5'/A^) (lower panels in the same boxes). 

The previous formulae allow us to remove the depen- 
dence of diagnostics on Rso and S/N . Consider now a galaxy 
for which we know the diagnostics Di and the quantities 
Rso/FWHM and S/N. For any given morphological type 
T and for each diagnostic, the deviations of the actual values 
of Di with respect to those obtained from the above fitting 
procedure, normalized by the corresponding expected scat- 
ter, can be expressed as follows: 



A - ^ ar'"(T) \og{R8o/FWHMr log(5/iV)^ 



AA(T) 



m+n=0,l,2 



If we assume the distributions of the deviates to be nor- 



for a galaxy of morphological type T, can be expressed by: 
P,{T) = SNF^[AD^{T)] 

(5'A^F[x]=e"^^/V^2^ is the Standardized Normal Func- 
tion) and the Maximum Likelihood (ML) probability that 
the morphological type T is associated to the actual set of di- 
agnostics Di (i=l,...,21) can be obtained through the prod- 
uct of the individual probabilities: 

P(T)= Yl P,{T) 

We actually prefer to use the logarithmic form: 
log[P(T)] =Y,'^i X log[Pi(T)] 

2 = 1,21 

where Wi are binary weighting factors (0/1) which determine 
the inclusion of each diagnostic in the final adopted set of 



Di (see Section 3.3.1 and Figure |8|. 

We compute the above probability from T=— 6.5 to 
T=10.5 (step 0.1) and we assume the ML morphologi- 
cal estimator to be the value of T which minimizes the 



— log [P(T)]. Figure B4 shows some examples of the be- 
haviour of the ML probability as a function of the morpho- 
logical type for typical galaxies belonging to the six above 
defined broad morphological classes. The red lines in each 
left side figure mark the minimum values of — log [P(T)], 
i.e. the adopted MORPHOT-ML morphological types Tml, 
while the dashed lines indicate the corresponding visual es- 
timates Tv . Conventionally, we should obtain the confidence 
intervals of our ML estimates by the change in T necessary 
to decrease P{T) by AP=0.5 from its value at the maximum 
(in logarithm: AlogP ^0.3). In practics, numerical simula- 
tions of the MORPHOT-ML technique and visual inspection 
of several real cases suggest that a value of AlogP ^0.2 
provides more realistic boundaries to the morphological es- 
timates. The green lines in Figure |B4| set the confidence 
intervals obtained in this way. Note in the figure the dis- 
continuity between cD and E galaxies already shown in Fig- 
ure |B2| and commented before. Note also the discrepancy 
between visual and ML estimates in the case of the galaxy 
in Abell 160, which has been visually classified early-spiral 
and automatically classified late-spiral (AT ^2.5). 



APPENDIX C: 
TECHNIQUE 



THE NEURAL NETWORK 



mal (a fair assumption, in our case; see Figure B3 ), the prob- 
ability that the actual value of the diagnostic Di is found 



The architecture and strategy adopted to produce a Neural 
Network estimator for the morphological classification are 
very similar to that described in jVanzella et al., ( ,2004| ). Here 
be briefly recall the method. 

CO. 1 Architecture 

We adopt the classical feed-forward multilayer perceptron 
Neural Network (MLP, Bishop fl 995 ) with four layers, each 
one made of 23, 20, 20 and 1 nodes. The first layer (named 
input layer) receives 23 input parameters for each galaxy 
(i.e. the diagnostics Di^ i=l,...,20 as described in Ap- 
pendix [A] and the global quantities £, \og[Rso / FW H M] and 
\og[S/N])^ while the single output node (output layer) pro- 
duces the morphological evaluation. The other two layers of 
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Spline interpolations for 
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Figure B2. Example of the fitting procedure performed in order to remove the dependence of morphological diagnostics on the size 
(Rso) and signal to noise ratio (S/N). The case of the diagnostic D13 is illustrated in the figure. Left panels: Natural Cubic Spline (NCS) 
fitting of the polynomial coefficients A^'^ and of the scatter ais as a functions of the morphological type (details in the text). Right 
panels: Residuals of the (polynomial+spline) fits illustrated in the left panels with respect to both log(i?8o/^^^^) and log(S'/A^). In 
particular, the upper panels inside each ATj tick line box show the residuals as a function of log(i?8o/^^^^) (units at the top of the 
figure), while the lower panels in each tick line box show the corresponding residuals as a function of log(»S'/A^) (units at the bottom of 
the figure). 



20 nodes each are called hidden layers. Each node of a layer 
is connected to all the nodes of the next layer. 



CO. 2 Training 

We use the supervised learning method, i.e. the NN is 
trained with examples. The training set is composed by the 
926 galaxies in the calibration sample, each one visually clas- 
sified (Tv). Each example presented to the NN is a pair of 
arrays, one contains the set of diagnostics (described in Ap- 
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Figure B4. Examples of the Maximum Likelihood minimization for WINGS galaxies of different 'broad' morphological classes. 



pendix [A| with 23 components and the second array is the algorithrrj^ modifies the strength of the connections between 
targeted morphology (single component), Ty . The learning 

^ Here we have used the hack-propagation algorithm with its gen- 
eralized delta rule version (see [Vanzella et al.|2004| ) 
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nodes (called weights) in order to force the association be- 
tween the input and the correct output. It translates to a 
minimization of a merit function as a function of the set of 
weights. 

Once a suitable set of weights is determined, it is frozen 
and stored for evaluation of new (never seen) galaxies. There 
are various techniques to identify the best set of weights, 
i.e. those that offer the so called best generalization power 
and that avoid the over- fitting problem. Here we adopt the 
method based on a committee of neural networks (30 MLPs 
have been used) as described in Vanzella et al.| ( [2004| . This 
method reduce the variance in the predictions maintaining 
a relatively small bias (Bishop 1995). Each member of the 
committee has been trained with different initial conditions 
(e.g. initial random distribution of weights, random sequence 
of examples presented to the NN, bootstrapping of exam- 
ples), that produce different histories of training and 30 dif- 
ferent set of weights. Therefore, for each galaxy, the commit- 
tee produces 30 estimations of the morphology (T^iv? ^=1- 
30). From this distribution we extract the median(mean) 
Tnn of the T^iv the central 68% interval (16 and 84 
percentiles). With this method we can associate a statistical 
uncertainty to the predicted value Tnn- 

It is clear that the outlined weighting of the connec- 
tions, inherent to this methodology, automatically provides 
the selection of the most significant diagnostics, thus mak- 
ing useless the preliminar (empirical) choice we performed 
in the case of the ML technique. 



